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g-breathers are exact time-periodic solutions of extended nonlinear systems continued from the 
normal modes of the corresponding linearized system. They are localized in the space of normal 
modes. The existence of these solutions in a weakly anharmonic atomic chain explained essential 
features of the Fermi-Pasta-Ulam (FPU) paradox. We study q-breathers in one- two- and three- 
dimensional discrete nonlinear Schodinger (DNLS) lattices — theoretical playgrounds for light prop- 
agation in nonlinear optical waveguide networks, and the dynamics of cold atoms in optical lattices. 
We prove the existence of these solutions for weak nonlinearity. We find that the localization of 
g-breathers is controlled by a single parameter which depends on the norm density, nonlinearity 
strength and seed wave vector. At a critical value of that parameter g-breathers delocalize via reso- 
nances, signaling a breakdown of the normal mode picture and a transition into strong mode-mode 
interaction regime. In particular this breakdown takes place at one of the edges of the normal 
mode spectrum, and in a singular way also in the center of that spectrum. A stability analysis of 
g-breathers supplements these findings. For three-dimensional lattices, we find g-breather vortices, 
which violate time reversal symmetry and generate a vortex ring fiow of energy in normal mode 
space. 

PACS numbers: 63.20.Pw, 63.20.Ry, 05.45.-a 



I. INTRODUCTION 
A. Background 

\ In 1955, Fermi, Pasta and Ulam (FPU) published their seminal report on the absence of 
I thermalization in arrays of particles connected by weakly nonlinear springs In particu- 
lar they observed that energy, initially seeded in a low-frequency normal mode of the linear 
problem with a frequency Uq and corresponding normal mode number g, stayed almost 
completely locked in a few neighbor modes in frequency space, instead of being distributed 
quickly among all modes of the system. The latter expectation was due to the fact, that 
nonlinearity does induce a long-range network of interactions among the normal modes. 
/^From the present perspective, the FPU problem appears to consist of the following parts: 
(i) for certain parameter ranges (energy, system size, nonlinearity strength, seed mode num- 
ber), excitations appear to stay exponentially localized in g-space of the normal modes for 
long but finite times [H ; (ii) this intermediate localized state is reached on a fast time scale 
Ti, and equipartition is reached on a second time scale T2, which can be many orders of 
magnitude larger than ri [sj; (iii) tuning the control parameters may lead to a drastical 
shortening of T2, until both time scales merge, and the metastable regime of localization is 
replaced by a fast relaxation towards thermal equilibrium |3|, which is related to the non- 
linear resonance overlap estimate by Izrailev and Chirikov FPU happened to compute 
cases with values of T2 which were inaccessible within their computation time. Among many 
other intriguing details of the evolution of the intermediate localized state, we mention the 
possibility of having an almost regular dynamics of the few strongly excited modes, leading 
to the familiar phenomenon of beating. That beating manifests as a recurrence of almost all 
energy into the originally excited mode (similar to the beating among two weakly coupled 
harmonic oscillators), and was explained by Zabusky and Kruskal with the help of real space 
solitons of the Kordeweg-de-Vries (KdV) equation, for a particular case of long wavelength 
seed modes, and periodic boundary conditions j^. The regular dynamics may be replaced 
by a weak chaotic one upon crossing yet another threshold in control parameters |6|. Re- 
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markably, that weak chaos is still confined to the core of the localized excitation, leaving 
the exponential localization almosts unchanged, and perhaps influencing most strongly the 
parameter dependence of T2. The interested reader may also consult Ref. 

B. Motivation 

Single-site excitations in translationally invariant lattices of interacting anharmonic os- 
cillators are known to show a similar behaviour of trapping the excitation on a few lattice 
sites around the originally excited one [sj. As conjectured almost 40 years ago [9|, ex- 
act time-periodic and spatially localized orbits - discrete breathers (also intrinsic localized 
modes, discrete solitons) persist in such lattices, and the dynamics on such an orbit, as well 
as the dynamics of the nearby phase space flow, account for many observations, and have 
found their applications in many different areas of physics [10]. Recently, it was shown, 
that the regime of localization in normal mode space for the FPU problem can be equally 
explained by obtaining q-breathers (QB) - time-periodic and modal-space-localized orbits 
- which persist in the FPU model for nonzero nonlinearity [ll|. The originally computed 
FPU trajectory stays close to such g-breather solutions for long times, and therefore many 
features of its short- and medium-time dynamics have been shown to be captured by the 
g-breather solution and small phase-space fluctuations around it. Delocalization thresh- 
olds for g-breathers are related to resonances, and corresponding overlap criteria The 
method of constructing g-breather solutions was generalized to two- and three-dimensional 
FPU lattices l3] . Most importantly, a scaling theory was developed, which allowed to con- 



struct g-breathers for arbitrarily large system sizes, and to obtain analytical estimates on 
the degree of localization of these solutions for macroscopic system size [l3|. Expectations 
were formulated, that the existence of g-breathers should be generic to many nonlinear 
spatially extended systems. 

C. Aim 

The spectrum of normal mode frequencies of the linear part of the FPU model is acoustic, 
contains zero, reflecting the fact, that the model conserves total mechanical momentum. 
In addition it is bounded by a lattice-induced cutoff, the analogue of the Debye cutoff in 
solid state physics. We may think of two pathways of extending the above discussed results 
to other system classes. First, we could consider a spatially continuous system, instead of 
a lattice. However, the initial FPU studies showed conflnement to a few long wavelength 
modes, and the results of Zabusky and Kruskal [s!], conflrm, that indeed in the spatially 
continuous system (KdV), trajectories similar to the FPU trajectory persist. 

As discussed above, the persistence, localization, and delocalization of g-breathers is due 
to a proper avoiding of, or harvesting on, resonances between normal modes. It therefore 
appears to be of substantial interest to extend the previous results to systems with a 
qualitatively different normal mode spectrum. Such a spectrum is one which does not 
contain zero, and is called optical. It may correspond to the excitation of certain degrees 
of freedom in a complex system, or to an FPU type model which is deposited on some 
substrate, or otherwise exposed to some external flelds. Together with the acoustic band, 
these two band structures, when combined, describe almost any type of normal mode 
spectrum in a linear system with spatially periodic modulated characteristics. 

Adding weak nonlinearities to such a system, taking local action angle representations, 
and performing various types of multiple scale analysis, such models are mapped on discrete 
nonlinear Schrodinger models (DNLS) These models enjoy a gauge invariance, and a 
conservation of the sum of the local actions (norms), i.e. some global norm. A particular 
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version of these equations is known as the discrete Gross-Pitaevsky equation, and is derived 
on mean fie ld g rounds for Bose-Einstein condensates of uhracold bosonic atoms in optical 



lattices [15|, Il6|. The norm is simply the conserved number of atoms in this case, and the 
nonlinearity derives from the atom-atom interaction. The propagation of light in (spatially 
modulated) optical media is another research area where DNLS models are used In 
that case the norm conservation derives from the conservation of electromagnetic wave en- 
ergy along the propagation distance (in the assumed absence of dissipative mechanisms). 
But DNLS models served equally well in many other areas of physics, where weak interac- 
tions start to play a role, e.g. in the theory of polaron formation due to electron-phonon 
interaction in solids. 

Releasing the norm conservation (e.g. by allowing the number of atoms in a condensate to 
fluctuate) will reduce the symmetries of the corresponding model, but most importantly, it 
will lead to possible new resonances of higher harmonics. We will discuss these limitations, 
and possible effects of releasing these limitations on g-breathers, in the discussion section. 

A particular consequence of norm conservation is a corresponding symmetry of the normal 
mode spectrum, which relates both (upper and lower) band edges, and makes the band 
center a symmetry point as well. 

The paper is structured as follows. Section [III is introducing the model and the main 
equations of motion. Section IIIII gives an existence proof for g-breathers. In section [IV] 
we derive analytical expressions for the QB profiles using perturbation theory. These are 
compared with the numerical results for d = 1 in section |Vl which also contains a stability 
analysis of the obtained periodic orbits. Section IVII gives results on periodic boundary 
conditions (as opposed to fixed boundaries). We generalize to o? = 2, 3 in section IVnl and 
discuss all results in section IVIIII 



II. MODEL 



We consider a discrete nonlinear Schrodinger (DNLS) equation on a (i-dimensional hy- 
percubic lattice of linear size N: 

meD{n) 

Here ^ is a complex scalar which may describe e.g. the probability amplitude of an atomic 
cloud on an optical lattice site [16|, or relates to the amplitudes of a propagating electromag- 
netic wave in an optical waveguide [13] • The lattice vectors n, m have integer components, 
and D{n) is the set of nearest neighbors for the lattice site n. If not noted otherwise, 
we consider fixed boundary conditions: ipn = if = or ri; = -|- 1 for any of the 
components of n. Equation ([1]) is derived from the Hamiltonian 

^ = e( E ^-c+fi^nr) (2) 

n \meD(n) j 

using the equations of motion z?/^„ = dH/dip^. In addition to energy, equation ([T]) conserves 
the norm -B = X] iV'nP- Note that the change of the nonlinearity parameter in ([T]) is 

n 

strictly equivalent to changing the norm B. Here we will keep the norm B (alternatively 
the norm density) fixed, and vary /x. 
We perform a canonical transformation to the reciprocal space of normal modes (g-space 
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of size A^'^) with new variables Qq{t) = Qq^...q^{t) 

^ r) \d/2 N d , s 

= (]vTt) E 0»Wn-(^)- (3) 

Together with ([T]) we obtain the following equations of motion in normal mode space: 

^Qq = —^qQq + /jy -j^w Cq^p^r,sQ pQ rQ (4) 



p,r,s 



where = —2 ^ cos -^^^ are the normal mode frequencies for the linearized system ([T]) 

i=l 

with /i = 0. Nonlinearity introduces a network of interactions among the normal mode 
oscillators with the following coupling coefficients: 

d 

Cq,p,r,s = ri ^qi,Pi,ri,Sii 
1=1 

^q^,P^,r^,s, = ( ~ ^ ) ('^(- + (- 1)' ^ + (- 1)'" s„g, + ^ ^ 

k,l,Tn=0 

+'^(-l)fcp,+(-l)'n+(-l)'"s„g,±(2Ar+2) ) • 

III. PROOF OF EXISTENCE OF t-REVERSIBLE g-BREATHER SOLUTIONS FOR WEAK 

NONLINEARITY 

We look for exact time-periodic solutions, which are stationary solutions of the DNLS 
equation ([1]), and which are localized in normal mode space: ipn{t) = (j)nexp{iQt) with 
frequency Q and time- independent amplitudes (pn- In the space of normal modes these 
stationary solutions have the form Qq(t) = Aqexp{iQt) , where the amplitudes of modes Aq 
are time- independent and related to the real-space amplitudes by the transformation ([3]). 
At a given norm B they satisfy a system of algebraic equations: 

— ^lAq + LOqAq — (^j^^-^^d XI ^ q,p,r,s^p^^^s = 0, (ll^..^d = '^1^1 

P,r,s (c\ 

E\Aq\'-B=0, 
Q 

We are focusing here (and throughout almost all of the paper) on t-reversible periodic orbits. 
Therefore we may consider all Aq to be real numbers. In this case system ([6]) contains + 1 
equations for + 1 variables. This system can be condensed into an equation for a vector 
function: 

F{X;fi,B) = (7) 

with X = {. . . , Aq, . . . , Q}. The components of F are the left hand sides of ([6]), while /i, 
B are parameters. 

For n = the normal modes in (jl]) are decoupled and each oscillator conserves its norm 
in time: Bq{t) = A"^. Let us consider the excitation of only one of the oscillators with the 
seed mode number Qq-. Bq = Bq^6qq^. The excited normal mode is a time-periodic solution 
of (jl]), and is localized in q-space. According to the Implicit Function Theorem [l8|, the 
corresponding solution of can be continued into the nonlinear case (/i 7^ 0), if the Jacoby 
matrix of the linear solution (5i^/(9X)^=o,B is invertible, i.e. ||c?i^/9X||^=o,B 7^ 0. The 
Jacobian 

=i-ir^'2AlU(-^'^-^J- (8) 

M=0 q^q^ 
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Invertibility of the Jacoby matrix requires a non-degenerated spectrum of normal mode 
frequencies. Therefore, the continuation of a hnear mode with the seed mode number Qq 
will be possible if ujq ^ cUq^, for all q ^ Qq. That condition holds for the case d = 1, and 
thus g-breathers exist at least for suitably small values of n there. For higher dimensions, 
degeneracies of normal mode frequencies appear. These degeneracies are not an obstacle 
for numerical continuation of g-breathers, but a formal persistence proof has to deal with 
them accordingly. An alog ous results for two- and three-dimensional /3-FPU lattices have 



been obtained in Ref. 12 



IV. PERTURBATION THEORY FOR g-BREATHER PROFILES 

To analyze the localization properties of g-breathers in g-space we use a perturbation 



)rop e 

theory approach similar to [IjJ, 12]. We consider the general case of a (i- dimensional DNLS 
lattice ([1]). Taking the solution for a linear normal mode with number Qq = (gi,o, • • • , Qdfi) 
as a zero-order approximation, an asymptotic expansion of the solution to the first 
equations of 1^ in powers of the small parameter a = ix2^^~'^^ /{N + 1)'^ is implemented. 
Note, that in the same way as it was done in 11], T2], we fix the amplitude of the seed mode 



Aq^. Later on we will apply the norm conservation law (the last equation of ([S])) to express 
the norm Bq^ = \Aq^^\'^ via the total norm B for the case of d = 1. Analytical estimations 
presented below describe amplitudes of modes located along the directions of the lattice 
axes starting from the mode Qq. Mode amplitudes of g-breathers have the slowest decay 
along these directions. Studying g-breather localization along a chosen dimension i, for the 
sake of compactness we will use a scalar mode number to denote the i-th component of q, 
assuming all other components be the same as in ^q: qj^i = qj Q. 

A. Close to the band edges 

Let us start with g-breathers localized in the low-frequency mode domain (g^ q << A^, 
i = l,d). According to the selection rules, if the seed mode number gjo is even (odd), 
then only even (odd) modes are excited along i-th dimension. The n-th order of the 
asymptotic expansion is the leading one for the mode gj_„ = (2n -|- l)gj^o- When gj reaches 
the upper band edge in the i-th direction, a reflection at the edge in g-space takes place if 
A^ + 1 is not divisible by g^^o (cf- Fig. [T]i). If A^ -|- 1 is divisible by gj_o then only modes 
qi^n = (2n -|- l)gi,o < A^ are excited. In the analytical estimates we assume a large enough 
lattice size. Then the effect of band edge reflections is appearing in higher orders of the 
perturbation, which will not be considered. In this case the approximate solution is 



^(2n+i).,o = (-sign(/i))" ^\/X']'j Aq^^,t = l,d,n = uq^^ - a A^^^^ + 0{a') . (9) 
The corresponding exponential decay of mode norms is 



^(2n+l)g,,o - (Ad ) Bq^, y (A^ ) - - 2 ' « " 1, « • (10) 

2^ '^TC^{qifi} 2^ [kifl) 

A^*'' (0 < A^*^ < 1) characterize the exponential decay along the i-th dimension. Here fcj^o = 
< TT is the i-th component of the seed wave vector ko, and fofeg = Bk^^/{N + lY 
is the norm density of the seed mode. When using intensive quantities - wave number fcj q 
and norm density 6fco - only, A does not depend on the system size. 

Equations (jl]) are invariant under the symmetry operation /x —fi, t — t, g^ q 
A^ -|- 1 — qifi for alH = 1, (i, which changes the sign of nonlinearity, and maps modes from 
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one band edge to the other. The replacement ^ N + 1 — for all i = l,d in ([6]) 
is equivalent to substitutions fi — > — and Q — > — fi. Using this symmetry, we can easily 
apply the above results to the case of g-breathers locahzed near the upper band edge, by 
counting mode indices from the upper edge: = iV + 1 — g^. We neglect reflections from 
the lower band edge, such that only modes with numbers = (2n + l)qifl {n is integer, 
qifl « N) are assumed to be excited. It follows 



/l;v+i-,~,„ = (sign(/i))" J K (11) 

where n is an integer and 

n-J^. _ \^.\Al^{N+lf-^ _ iMlbfeQ . _ j-^ (12) 



The analytical expression for the frequency of the g-breather solution turns out to be the 
same as for the case of small seed mode numbers. 

B. Close to the band center 

We implement the perturbation theory approach for seed modes close to the band center. 
Let us first consider the case of odd A^. We introduce the new index p\ which is the 
number of a mode counted from the middle of the spectrum along the i-th dimension: 
Pi = Qi — {N + l)/2. We choose the seed mode with \pifl\ « (iV + l)/2. In the n-th order 
of perturbation theory the newly excited mode along the i-th dimension has the number 
Pi,n = (-l)"(2n + l)pifi, and 

= (-A?)M,^,o , = signfe,o)(-l)"+^(A?)"+^/^A,,,o , (13) 

where i = l,d, A^^ > 0, n = 1, 2, 3 and 



(14) 



24"'^vr|g,,o - (A^ + 1)/2| 2^~''\h,o - vr/2| ' 

When using intensive quantities, A again does not depend on the system size. 

For the case of even N we use pi = qi — N/2 and assume that the seed mode index pj^o > 0, 
i = l,d. The set of consecutively perturbed modes becomes more complicated: 

Plo ^ (-3pi,o + 2) (5pi,o-2) ^ (-7pi,o + 4) (9p,,o-4) ^ (-11^-0 + 6) ^ ... . (15) 

For an even number of perturbation steps the new excited mode along the i-th dimension 
has index pi^2n = (4n + l)pifi — 2n, and for an odd number of steps it has index Pi,2n+i = 
— (4n + 3)pifi + 2n + 2, n = 0, 1, 2, . . .. The amplitudes satisfy (fT3|) . but with 

(16) 

For large enough equation ( JT6ll approaches the expression (IHl), therefore we will use (JT^ 
for large only. 

V. g-BREATHERS FOR d = 1: RESULTS 
A. Close to the band edge 











l^fco 


2^-'^vr|gi,o - iV/2| 2^"'^\ 


\ki,o - vr 


/2 + 7r/(2Ar + 2)| 



The key property of g-breathers is that they are localized in the space of linear normal 
modes. Note that some g-breather solutions may be compact in g-space and contain only 
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FIG. 1: Distributions of mode norms for QBs in the one-dimensional DNLS model with parameters: B = 1, fi — 0.1 
for (a), (b), (c); (a) QBs with low- frequency seed mode qo = 1, (b) QBs with high-frequency seed mode qo = N, (c) 
QBs with seed mode near the middle of the spectrum qo = H- (A'^-l- 1)/2 (c). Dashed lines in (a), (b), (c) correspond 
to analytical estimations, (d) Observation of multiple reflections at the boundaries of the g-space, A'^ = 31, B — 1, 
fj, = 2, qo = 3. Dashed lines in (d) are guidelines for the eye. 
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FIG. 2: The slope S as a function of the seed wave number ko for d — 1 with fixed norm density b — 1/1025, 
fi — ±0.5, ±1, ±2 (from bottom to top). Dashed lines correspond to the analytical estimates. Different symbols: 
slope from the numerical calculation of QB (squares and circles represent the results for A'' — 512). Solid lines guide 
the eye. 



Let US consider the one- dimensional case. We compute g-breathers as the stationary 
solutions of the nonlinear equations ([6]) using the single-mode solution for /i = as an 
initial approximation. Fig. [T^,b shows the distribution of mode norms for g-breathers in 
the space of wave numbers k = 7Tq/{N + 1) for different chain sizes and different seed wave 
numbers, located near the lower {ko « n, Fig. [T^) and upper (vr — feg << tt, Fig. [T)d) 
edges of the linear mode spectrum. 

We find exponential localization of g-breathers, with a localization length which depends 
strongly on the chosen parameters. The obtained analytical estimations for g-breathers 
localized at the lower band edge fl9|). f|T0|) respectively upper band edge f|TTl) . f|T2|) . are in 
quantitative agreement with numerical results (see dashed lines in Fig. [T^,b). 

The exponential decay in ffTOl) is depending on the seed mode norm density. Many appli- 
cations (e.g. cold atoms in a condensate) rather fix the total norm, or total norm density. 
For the obtained g-breather solutions with qq << N, the relation between these quantities 
can be estimated by the sum of infinite geometric series: 

oo „ 

5^J]A"5,„ = ^, (17) 

n=0 

where A = X^^\ From ([17]) it follows that = (1 - A)6, where b = B/{N + 1) is the 
total norm density. Substituting the expression for into (ITOl) and solving the equation 
for -\/A we obtain: 



(18) 
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The same dependence of on ko was obtained for low frequency g-breathers in the /3-FPU 



model [13|. The exponential decay of mode norms in the space of wave numbers, can be 



now written for g-breathers localized in low-frequency modes: 



\nbk= {^-l)\nVX + \nbk,. (19) 



ko 

To characterize the degree of localization in /c-space, we use the slope of the profile of mode 



norms in log- normal plots ffT^ - S [13|, where the absolute value of S is equivalent to the 
inverse localization length ^: 

S=^\n^/X, \S\^r'- (20) 

Substituting the expression for (UHl) into ([20]) we obtain: 

5 = -1 In v^l + zyi - , z = ko/u . (21) 

Therefore the slope (inverse localization length) is a function of the rescaled wavenumber z. 
It therefore parametrically depends on just one effective nonlinearity parameter u, which 
is given by the product of the total norm density and the absolute value of nonlinearity 
strength. S vanishes for z ^ 0, and it has an extremum min(S') ~ — 0.7432/z/ at z^in ~ 
2.577. The wave number ko = kmin = -Zmm'^ corresponds to the strongest locahzation of a 
g-breather with fixed effective nonlinear parameter v. With increasing v the localization 
length increases. Most importantly, the localization length diverges for small k^ <^ kmin 
since l^l ~ z/{2v) in that case. This delocalization is due to resonances with nearby normal 
modes close to the band edge. Note, that our analytical estimations do not depend on the 
sign of the nonlinearity parameter /x. 

In Fig. [T^,b, due to the small sizes of the chain, all values of k^ are greater than kmin-, 
so a monotonous dependence of the slope on k^ is observed. In Fig. [2] we plot theoretical 
and numerically obtained dependencies S{kQ) for different values of the nonlinear param- 
eter /i and different system sizes N (we use large enough to resolve the theoretically 
predicted extremum of S). In all cases the numerical results show, that the slopes indeed 
are characterized by intensive quantities only, and the above derived scaling laws hold. 

For positive values of and seed wave numbers close to the lower band edge, the extremum 
in S is reproduced in the numerical data, though the numerical curves deviate from theo- 
retical curves for small /cq- We also find that k^m increases and \S{kmin)\ decreases with 
increasing /i as it is predicted by our analytical results (increase of norm density h gives the 
same effect). 

For negative values of /i and seed wave numbers close to the lower band edge, we do 
not observe an extremum for S. In the region of small ko the numerically obtained curves 
for positive and negative values of differ, while the theoretically predicted slopes do not 
depend on the sign of /i. Fig. [3] shows the mode norm profiles of g-breather solutions with 
small ko obtained for positive and negative values of /x. The reason for the discrepancy be- 
tween theoretical prediction and numerical results for negative values of /x must be strong 
contributions from higher order terms in the perturbation expansion. The standard argu- 
ment is, that the perturbation theory is valid if the localization length is small, i.e. |A| ^ 1 
as well. When |A| becomes of the order of one, higher order terms in perturbation theory 
have to be taken into account, and it would be tempting to conclude that delocalization 
will take place. That should be true especially when all higher order terms in the series 
carry the same sign. That is the case for positive nonlinearity here, but for negative we 
obtain alternating signs of higher order terms. These alternating signs therefore effectively 
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FIG. 3: Distributions of mode norms for QBs in the one-dimensional DNLS model for positive and negative values 
of At; go = 2, iV = 1024. 
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FIG. 4: The master slope function Sm{z) (dashed line). Different symbols and eye-guiding solid lines correspond to 
the scaled numerical estimates of the slope presented in Fig. [5] (bottom line — A* < Oi top line — > 0). 



cancel most of the terms in the series, and are responsible for an increasing localization of 
g-breathers with negative n in the limit of small wavenumbers. 

The obtained results for g-breathers with seed wave numbers close to the lower band 
edge {ko « n) are valid for the case of ko close to the upper band edge (vr — /cq << vr) 
if we change fi —fi. Thus, there is a strong asymmetry in the localization properties of 
g-breather solutions with ko < 7t/2 and ko > 7r/2 for a fixed sign of nonlinearity. 

The numerical results in Fig. [2] show that slope values calculated for different system 
sizes lie on the same curves with corresponding /z even in the region of small ko, where 
higher order corrections to our analytical estimates have to be taken into account. This 



result is in agreement with the exact scaling of g-breather solutions described in [13|. We 



plot in Fig. mthe master slope function Sm{z) = uS, which depends on a single variable 



11 



z [13|. It implies that knowing this single master slope function is sufficient to predict the 
localization property of a g-breather at any seed wave number <^ vr, at any energy etc. 
Numerically obtained slopes presented in Fig. [2] are rescaled and plotted in Fig. HI We see, 
that all results corresponding to the same sign of condense on a single curve even (and 
especially) for small /cq, though these numerical results differ from the analytical estimation. 



B. Close to the center of the band 



Fig. [T}:; shows the distribution of mode norms for g-breathers in the space of wave numbers 
k = TTq/{N + 1) for different chain sizes and different seed wave numbers, located close to 
the center of the spectrum (|/co ~7r/2| << vr). The analytical estimation for the amplitudes 
of these g-breathers are in good quantitative agreement with the numerical results, for small 
enough parameters /i and 6, cf. Fig. [Tt. 

We express again the norm density of the seed mode via the total norm density: fo^p = 
(1 — A)6, (A = A^^'*). Substituting this expression into ffTll) we find: 

/y _ v/l+4i/V(fco-7r/2)2-l _ |^ /22^ 
y — 2u/\kQ-K/2\ ' ~ 8 • ^ ' 

The slope of the mode norm proffie in fc-space is given by 



S 



1^0 -7r/2 



In = — In f VlT^ - z\ , (23) 
2i'z V / 



where z = \kQ - 7r/2|/(2z/). For z « I: S = (-1 + z^/Q + 0{z^))/{2iy), for z » 1: 
S = ln(l/(22) + 0{l/z^))/{2uz). In the limit z ^ 0, the slope S -l/(2i/). Thus, the 
strongest localization of a g-breather with fixed effective nonlinear parameter u should be 
obtained for wave numbers ko ~ 7r/2. The increase of the effective nonlinearity parameter 
u ~ \fi\b leads to a weaker localization of g-breathers in fc-space, since the absolute value 
of the slope S decreases. Note, that there is a special point for odd N: ko = vr/2, go = 
(A^ + l)/2. For this seed mode, the g-breather is compact in fc-space (l9l |. 

In Fig. Elthe theoretical and numerically obtained dependencies S{ko) for different values 
of fi and different system sizes are plotted. For small values of /i we observe good 
agreement between analytical and numerical results. But the increase of the nonlinearity 
leads to a deviation between the theoretical and numerical curves in the region of ko close 
to 7r/2. These corrections, as it was for the case of g-breathers localized near the band 
edges, depend on the location of ko (fco > 7r/2 or fco < vr/2) and on the sign of nonlinearity. 
Therefore the curves of S{ko) for strong nonlinearity {fx = 30) in Fig. [5] are non-symmetric 
around the point ko = n /2: in contrast to ko < 7r/2, for ko > 7r/2 a local minimum of 5* 
is observed. Still, the predicted scaling properties of g-breathers remain correct even for 
strong nonlinearity: the values of S, computed for different system sizes A^, lie on the same 
curves for fixed /x. 



C. StabiHty of g-breathers 



We analyze the linear stability of g-breathers as stationary solutions of DNLS model con- 
sidering the evolution of small perturbations e„ in the rotating frame of the periodic solution 
(2^ : ipn{t) = (0^ + £^Ti(^)) exp(2r2t), where 0° are the non-perturbed time-independent am- 
plitudes. En = an + iPn- A pcriodic orbit is stable when all perturbations do not grow in 
time. Solving the linearized equations for the perturbation, that condition translates into 
the request, that all eigenvalues Sm {jn = 1, 2A^'^) of the linearized equations must be purely 
imaginary. Otherwise the orbit is unstable. In Fig. [6] we plot the numerical outcome of the 
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FIG. 5: The slope S as a function of the seed wave number ko for d — 1 with fixed norm density b = 1/1025, 
fi = 0.5,5,30 (from bottom to top). Dashed lines correspond to the analytical estimates. Different symbols 
correspond to the estimates of the slopes from numerical calculations of QBs. 



stability analysis. We smoothly continue g-breather solutions for each seed mode number 
go by increasing the nonlinearity parameter fi and check the stability. If the maximum 
absolute value of the real parts of all eigenvalues is smaller than 10~^, the g-breather is 
considered as stable (solid line), otherwise it is unstable (dashed line), crosses mark the 
change of stability. For small values of nonlinearity all g-breathers are stable. Qualitatively 
different threshold values and dependencies on go for g-breathers with seed modes from 
different parts of the spectrum gg < N/2 and go > N/2 are observed. This is in agreement 
with the results presented in [2l|, where stability properties of nonlinear standing waves 
are studied. Note, that due to the symmetry of the equations, the stability properties of 
a g-breather with seed mode go for some negative value of the nonlinearity parameter /x is 
the same as the stability property of the g-breather with seed mode go = + 1 — g for the 
nonlinearity parameter Jl = —fi. 

Let us discuss the possible link between linear stability and localization. If a g-breather 
becomes delocalized, that happens because of resonances between different mode frequen- 
cies. Therefore we can expect, that the same resonances will also drive the state unstable. 
Indeed, these correlations can be clearly observed from the numerical data. However, if a 
g-breather is well localized, it does not follow that it will be stable as well, since instability 
can arise due to resonant interaction of modes in the breather core alone. 
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FIG. 6: Domain of stability of g-breathers for d = 1. Solid lines - stable, dashed - unstable, crosses - switch from 
stable to unstable or vice versa. A'' = 64, B = 2. 



VI. PERIODIC BOUNDARY CONDITIONS 



In the case of periodic boundary conditions, we have used the following transformations 
between real space and the reciprocal space of normal modes (for even A^): 

Mt) = j^ E Q.W n exp (^^^^) . (24) 

qu...,qd=-N/2+l i=l ^ ^ 

The DNLS model ([1]) with periodic boundary conditions has exact solutions for nonlinear 
traveling waves, which can be written, for instance in case d = 1, as: tpn{t) = 0oexp(if2t — 
ikon), where Q = —2 cos ko — ficpl, ko = 27iqo/N, go G [—N/2, N/2]. These types of solutions 
can be considered as compact g-breathers which contain only one mode go- Traveling modes 
are also not invariant under time reversal. The continuation of a linear standing wave, 
consisting of two traveling waves with the same norms and wave numbers: fco and — /cq, into 
the nonlinear regime leads to a time- reversible g-breather solution (see Fig. [7]), which is 
not compact, and its localization properties are similar to the properties of g-breathers in 
the case of fixed boundary conditions. Here we present the result for decay of mode norms 
A^*'' for the case {h^l << it, which differs from (fTOl) by a prefactor: 

where fefcg = B^^^/N'^, ki^ = 27iqifi/N. Fig. [7| illustrates good quantitative agreement 
between analytical and numerical results obtained for small enough values of norm and 
nonlinearity. 
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FIG. 7: Time-reversible g-breather in the one-dimensional DNLS model with periodic boundary conditions, con- 
tinued from the linear standing wave consisting of two traveling waves with go = 2 and go ~ —2, A*' = 64, B — 1, 
fj, — 0.1. Dashed lines represent analytical estimations. 



VII. g-BREATHERS IN TWO- AND THREE-DIMENSIONAL LATTICES 

For two- and three-dimensional symmetric DNLS lattices ([1]) with fixed boundary condi- 
tions, only linear modes with mode numbers q on the main diagonal have non-degenerate 
frequencies. Using the Implicit Function Theorem [l8|, it follows that these modes are 
continued into the nonlinear regime. However, we also successfully performed numerical 
continuations of g-breathers with seed mode numbers off the main diagonal as it was done 
in tlH for the FPU model. For the two-dimensional DNLS model the g-breather, continued 
from the single linear mode Qq = (2,3), is presented in Fig. [8] in real space (a) and in 
g-space (b). The g-breather with the seed mode Qq = (3, 2) exists as well and has the same 
frequency. The slowest decay of the mode norms happens to be along the direction of the 
main axes. The decay is exponential and in good agreement with the analytical estimation 
( ITOj) for d = 2. In addition to such asymmetric single mode g-breathers it is possible to 
construct multi-mode g-breathers continued from a pair of degenerate linear normal modes 
(go,Po) and (j5o,go) (go 7^ Po) with the same norms in both modes and in-phase (Fig. [St) 
or antiphase (Fig. [Hli) oscillations. It is important to note that the problem of degenerate 
frequencies is avoided for these solutions. Indeed, system (jl]) has two invariant manifolds 
Qq^,q2 = ±<5g2,iji- Lookiug for a solution on a manifold, the number of independent variables 
of state is reduced from A^^ to the dimensionality of the manifold, which equals (A^^ + A^) /2 
for the symmetric manifold and [N^ — N)/2 for the antisymmetric one. The reduced system 
of equations contains only modes with non-degenerate frequencies. 

For d = 3 we have also verified, that the analytical estimations of mode norm decay flTU]) . 
f[T^ agree well with the results of numerical calculations for single-mode g-breathers. 

In addition to various time-reversible g-breather solutions, which are constructed in the 
same way as for = 2, the three-dimensional DNLS model allows also for non-time- 
reversible ("vortex") multi-mode solutions. Let us consider g-breather solutions on an 
invariant manifold of the system (jlj): Qq^,q2,q-j, = exp(27ri/3)(5q3,gi,g2 = 6xp(47ri/3)(5g2, 93,91 
(note, that this manifold has a counterpart with the opposite sign of the phase shifts). On 
the manifold the number of variables of state is reduced to (A^'^ — A^)/3. We have constructed 
numerically a vortex g-breather solution continued from a degenerate triplet of seed modes 
which have the same norm and a relative phase shift 27r/3: Qgo.qcpo — 6xp(27r2/3)QpQ^g(,^g(, = 
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FIG. 8: Different g-breather modes for d — 2: N = 32, B = 1, fi = 0.5. Mode norm magnitude is plotted in 
color code on logarithmic scale (except for (a), where a linear scale is used), (a), (b) — asymmetric single- mode 
g-breather with seed mode (qq = (3,2)) in real space (a) and in q-space (b); (c), (d) — symmetric multi-mode 
q-breather with in-phase (c) and antiphase (d) seed mode pair q = (3, 1), (1,3). Dashed line in (b) guides the eye 
along the main diagonal. 

exp(47ri/3)(5qo,po,go; 1o Po- The frequency of such a triplet becomes non-degenerate in the 
reduced system on the manifold. The energy flows in a vortex-like manner in g-space for 
such excitations (FigjQj). 

VIII. DISCUSSION 

Comparing the aims we formulated in the introduction with the above results, we can 
conclude, that indeed the g-breather concept turns out to be generic, and time-periodic or- 
bits, which are localized in normal mode space, are probably as generic as discrete breathers 
(although perhaps in a different way). Along with the study of the details of g-breathers 
properties in DNLS models, we found some particularities, which were not (yet) obtained 
in acoustic (FPU) systems. Let us discuss some of these and other findings from above. 
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FIG. 9: A vortex q-breather state for d — 3: N — 8, B — 1, jj, — 0.1, at some fixed time. Mode norm magnitude 
is encoded in sphere sizes on logarithmic scale, while the sphere color denotes the phase of the mode. Seed mode 
triplet: q = (3, 1, 1), (1, 3, 1), (1, 1, 3). 



A. Scales, delocalization thresholds, and healing length of a BEC 

Let us consider a finite DNLS chain with sites. Let us fix the total norm density b 
and the nonhnearity fi. Therefore we fix the effective nonhnearity parameter z/^ = |/i|6/16 
f|T8|) . If the system size was large enough, we will resolve the extremum in S, and therefore 
the QB with seed wave number fc^m ~ l.'bHv is the strongest localized one. It evidently 
sets an inverse length scale v. This length scale is known in the GP equation for a BEC in 
a trap with (similar) fixed boundary conditions. One has a condensate, whose amplitude 
vanishes at some boundary, yet getting back to some mean value away from the boundary 
- exactly at the healing length = (47ran)~^/^, where n is the condensate density, and a is 
the scattering length which is proportional to the atom-atom interaction, and therefore to 



the nonhnearity strength /z within the mean field GP equation [16|. Therefore, the inverse 
healing length corresponds to the wave number scale, on which the most strongly localized 
QB is observed. 

We now increase the system size further, and compute the localization length (or respec- 
tively its negative inverse - the slope S") of the longest wavelength mode. With increasing 
the grid of allowed /c-values becomes denser, and at some critical Nciy) the localization 
length will reach the finite size of the normal mode space, and the g-breather delocalizes. 
With a little algebra it follows 

iV, ^ 1^ , fcf) ^ ^ . (26) 
For even larger (and finally infinitely large) lattices the critical value k^^^ is marking a 
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border in /c-space: at the given norm density z/, all QBs with seed wave numbers k ^ /cq 
are localized, while we obtain delocalization for k < k^\ So there is a layer of delocalized 
QBs at the band edge, whose width grows according to (1261) with growing u. Modes 
launched inside this layer (with the given norm density) will quickly spread their energy 
among many other modes - they will quickly relax. Modes launched outside this layer will 
stay localized in normal mode space, at least for sufficiently long times, k^^ is therefore 
separating a layer of strongly interacting modes from weakly interacting ones. For large 
enough effective nonlinearity parameter i/ ~ 1 the whole wave number space is filled with 
strongly interacting normal modes, and the normal mode picture breaks down completely. 
As long as u is smaller, the value of k^'^^ sets a new length scale 2n/k^'^\ which is proportional 
to the squared healing length C,^. On that new length scale, the normal mode picture breaks 
down. 

B. Delocalization thresholds and modulational instability 

It is instructive to remember, that delocalization thresholds of QBs are related to reso- 
nances. Indeed, in the present case, the density of states at the band edge diverges in the 
limit of large system size, and many modes have almost the same frequencies. It is these 
small differences, which tend to zero as l/N"^, and which are responsible for the resonant 
mode-mode interaction. Notably the analysis of stability of band edge modes [23I] (note: 
for periodic boundary conditions) shows some interesting correlations. The analyzed band 
edge modes are compact in normal mode space, yet they undergo a tangent bifurcation at 
amplitudes, which are the smaller, the larger the system size. These instabilities appear 
however only for a certain sign of the nonlinearity, which exactly corresponds to the actual 
observed delocalization of QBs. Therefore we expect, that the (so far not studied) case of 
a FPU chain with negative quartic nonlinearity, which is known not to yield an instability 
for the band edge mode, will not show delocalization of QBs close to the (upper) band 
edge. It is furthermore instructive, that if a tangent bifurcation of a band edge mode takes 
place, the simulation of a lattice shows the onset of modulational instability, which leads 
to a collection of energy in smaller system volumes, and finally to the formation of discrete 
breathers - i.e. to localization in real space. Thus we may expect, that the resonant layer 
of strongly interacting modes may lead to the formation of localized states in real space, 
while the rest of the normal mode space is evolving in the regime of localization in normal 
mode space. Thus, we may expect to observe in an actual simulation localization both in 
normal mode space and in real space. 

C. Comparing analytical and numerical results 

While the delocalization at one band edge, as predicted by perturbation theory, is repro- 
duced by numerical data, that does not happen at the second band edge (both band edges 
change their places, when the nonlinearity inverts sign). The breakdown of perturbation 
theory follows from comparing the leading and next-to-leading order terms in the expansion. 
When terms are of the same order, we conclude, that the series will diverge, and therefore 
the QB solutions will delocalize. This is true for the case when all terms in the series have 
the same sign. However, when the terms have alternating sign, the above conclusion must 
not be correct. And indeed the numerical data show, that these sign alternations lead to an 
effective cancelation, and a final strong localization of a QB solution. A similar breakdown 
of perturbation theory happens for QB solutions localized at the band center for large non- 
linearity (or norm). The perturbation theory tells, that QB states are well localized from 
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both sides of the center. Numerical data however show, that this is true only on one side 
from the center, while the other side shows a tendency towards delocalization. 



D. Open questions 

Below we list some potentially interesting and important open questions. First, the 
influence of different boundary conditions has not been systematically studied. QBs are 
extended states in real space, therefore their spectrum, and the way they interact, may to 
some extend be sensitive to the choice of boundary conditions. 

Second, it remains completely open, what kind of excitation we catch by considering a 
vortex state in normal mode space. 

Finally, it would be interesting to see the changes in the QB properties, when the norm 
conservation is lifted. In general we expect, that this will lead to the generation of higher 
harmonics, and new types of resonances. In fact it is exactly resonances due to higher 
harmonics, which give the leading order contribution to the FPU problem. However, for 
sufficiently narrow optical bands, higher harmonics will be located outside the optical band. 
These resonances will be therefore important, when the optical band is wide, or when one 
considers a complex and broader band structure. 
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